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(N . Abstract 

o i 

I We explore the phase diagram of strongly interacting matter with 

' massless u and d quarks as a function of the strange quark mass m s and the 

Q\ ■ chemical potential \i for baryon number. Neglecting electromagnetism, we 

describe the different baryonic and quark matter phases at zero temperature. 
For quark matter, we support our model-independent arguments with a 
quantitative analysis of a model which uses a four-fermion interaction 
Q-i, abstracted from single-gluon exchange. For any finite m s , at sufficiently 

| large fi we find quark matter in a color-flavor locked state which leaves a 

global vector-like SU(2) co \ ov+ x j+ ji symmetry unbroken. As a consequence, 
chiral symmetry is always broken in sufficiently dense quark matter. As 
rS [ the density is reduced, for sufficiently large m s we observe a first order 

transition from the color-flavor locked phase to a color superconducting 
phase analogous to that in two flavor QCD. At this unlocking transition 
chiral symmetry is restored. For realistic values of m s our analysis indicates 
that chiral symmetry breaking may be present for all densities down to those 
characteristic of baryonic matter. This supports the idea that quark matter 
and baryonic matter may be continuously connected in nature. We map 
the gaps at the quark Fermi surfaces in the high density color-flavor locked 
phase onto gaps at the baryon Fermi surfaces at low densities. 



1 Introduction and Phase Diagram 

Strongly interacting matter at high baryon number density and low temperature is 
far less well understood than strongly interacting matter at high temperature and 
zero baryon number density. At high temperatures, the symmetry of the lowest free 
energy state is not in dispute, calculations using lattice gauge theory are bringing 
the equilibrium thermodynamics under reasonable quantitative control, and even 
non-equilibrium dynamical questions are being addressed. At high densities, we 
are still learning about the symmetries of the lowest free energy state. In addition 
to being of fundamental interest, an understanding of the symmetry properties 
of dense matter can be expected to inform our understanding of neutron star 
astrophysics and perhaps also heavy ion collisions which achieve high baryon 
densities without reaching very high temperatures. 

At high densities and low temperatures, the relevant degrees of freedom are 
those which involve quarks with momenta near the Fermi surface(s). The presence 
of an arbitrarily weak attraction between pairs of quarks results in the formation 
of a condensate of quark Cooper pairs. Pairs of quarks cannot be color singlets, 
and in QCD with two flavors of massless quarks the Cooper pairs form in the color 
3 channel.]]], 0, |3|, f|, The resulting condensate gives gaps to quarks with two 
of three colors and breaks the local color symmetry SU(3) co \ OI to an SU(2) color 
subgroup. The breaking of a gauge symmetry cannot be characterized by a gauge 
invariant local order parameter which vanishes on one side of a phase boundary. 
The superconducting phase can be characterized rigorously only by its global 
symmetries. In QCD with two flavors, the Cooper pairs are ud — du flavor singlets 
and, in particular, the global flavor SU{2) L x SU{2) R symmetry is left intact. 
There is an unbroken global symmetry which plays the role of baryon number 
symmetry, U(1)b- Thus, no global symmetries are broken and the only putative 
Goldstone bosons are those five which become the longitudinal parts of the five 
gluons which acquire masses. HQ 

In QCD with three flavors of massless quarks the Cooper pairs cannot be 
flavor singlets, and both color and flavor symmetries are necessarily broken. The 
symmetries of the phase which results have been analyzed in Ref. ||. The 
attractive channel favored by one-gluon exchange exhibits "color-flavor locking". 
It locks SU(3)l flavor rotations to SU(3) co \ OT , in the sense that the condensate 
is not symmetric under either alone, but is symmetric under the simultaneous 
5'L r (3) co i or+ z / rotations. The condensate also locks SU(3)r rotations to SU(3) co \ OT , 
and since color is a vector symmetry the chiral SU(3)l-r symmetry is broken. 
Thus, in quark matter with three massless quarks, the SU(3) co \ or x SU(3)l x 

1 There is also an unbroken gauged symmetry which plays the role of electromagnetism. Also, 
the third color quarks can condense, || as we discuss below, but the resulting gap is much smaller. 
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SU(3)r x U(1)b symmetry is broken down to the global diagonal <S77(3) co ior+v 
subgroup. There is also an unbroken gauged U(l) symmetry (under which all 
quarks have integer charges) which plays the role of electromagnetism. All nine 
quarks have a gap. All eight gluons get a mass. There are nine massless Nambu 
Goldstone excitations of the condensate of Cooper pairs which result from the 
breaking of the axial SU(3)a and baryon number U(1)b- We see that cold dense 
quark matter has rather different global symmetries for m s = than for m s = oo. 

A nonzero strange quark mass explicitly breaks the flavor SU(3)y symmetry. 
As a consequence, color-flavor locking with an unbroken global SU(3) co \ or+ v occurs 
only for m s = 0. Instead, for nonzero but sufficiently small strange quark 
mass we expect, and find, color-flavor locking which leaves a global SU(2) co \ or+v 
group unbroken. As m s is increased from zero to infinity, there has to be 
some value m^ nlock at which color and flavor rotations are unlocked, and the full 
SU (2) l x SU (2)r symmetry is restored. We argue on general grounds in Section ^ 
that this unlocking phase transition must be of first order. In subsequent Sections, 
we analyze this transition quantitatively in a model using a four-fermion interaction 
with quantum numbers abstracted from single-gluon exchange. 

From our analysis of the unlocking transition, we conclude that for realistic 
values of the strange quark mass chiral symmetry breaking may be present for 
densities all the way down to those characteristic of baryonic matter. This 
raises the possibility that quark matter and baryonic matter may be continuously 
connected in nature, as Schafer and Wilczek have conjectured for QCD with three 
massless quarks. We use our calculations of the properties of color-flavor locked 
superconducting strange quark matter to map the gaps due to pairing at the quark 
Fermi surfaces onto gaps due to pairing at the baryon Fermi surfaces in superfluid 
baryonic matter consisting of nucleons, A's, E's, and H's. (See Section^). 

We argue that color-flavor locking will always occur for sufficiently large 
chemical potential, for any nonzero, finite m s . We make this argument by first 
using our results as a guide to quark matter at moderate densities and then using 
them to normalize Son's model-independent analysis valid at very high densities. || 
As a consequence of color-flavor locking, chiral symmetry is spontaneously broken 
even at asymptotically high densities, in sharp contrast to the well established 
restoration of chiral symmetry at high temperature. 

Our goal in this paper is a consistent picture of the symmetries of cold dense 
quark matter with massless u and d quarks as a function of the strange quark mass 
m s and the chemical potential \i. We take /i to be the chemical potential for quark 
number, one third that for baryon number. m s is the current quark mass; we will 
refer to the //-dependent constituent quark mass as M s (/i). 

Figure 1 summarizes our conjecture for the zero temperature phase diagram of 
QCD as a function of m s and \i. (We work at zero temperature throughout this 
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Figure 1: Conjectured phase diagram for QCD with two massless quarks and 
a strange quark at zero temperature. The global symmetries of each phase are 
labelled. The regions of the phase diagram labelled 2SC, 2SC+S and CFL denote 
color superconducting quark matter phases. The Figure is described at length in 
the text. 



paper.) The rest of this section can be read as a description of this Figure. Lines in 
the diagram separate phases which differ in their global symmetries. This means 
that each line describes a distinction which can be associated with a local order 
parameter which vanishes on one side of the line. In each region of the diagram, we 
list the unbroken global symmetries of the corresponding phase. We characterize 
the phases using the SU{2) L x SU{2) R flavor rotations of the light quarks, and the 
U(l)s rotations of the strange quarks.Q 

In subsequent Sections of this paper, we explore those phases (labelled 2SC, 
2SC+s and CFL) which extend to high enough density that they can be described 
as superconducting quark matter. We analyze these phases, and the unlocking 
phase transition which separates them, quantitatively in a model in which quarks 

2 The U{1)b symmetry associated with baryon number is a combination of U(l)s, a U{\) 
subgroup of isospin, and the gauged C/(1)em of electromagnetism. Therefore, in our analysis 
of the global symmetries, once we have analyzed isospin and strangeness, considering baryon 
number adds nothing new. 
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interact via a four-quark interaction modelled on that induced by single-gluon 
exchange. Certainly, our model is not expected to be a valid description for QCD 
at nuclear matter density, where confinement plays an important role. That part 
of the diagram which describes the symmetry properties of different phases of 
baryonic matter is therefore not derived from our model. It is conjectural but 
plausible. 

In Figure 1 and throughout this paper, we neglect the small u and d 
current quark masses. The light quark masses have no substantial influence 
on the condensation of quark Cooper pairs. |5|, |9|] We also ignore the effects of 
electromagnetism throughout. We assume that wherever a baryon Fermi surface 
is present, baryons always pair at zero temperature. To simplify our analysis, we 
assume that baryons always pair in channels which preserve rotational invariance, 
breaking internal symmetries such as isospin if necessary. 

We now explain the features shown in Figure 1 by beginning at \x = (in 
vacuum) and describing the phase transitions which occur as /i is increased at 
constant m s . We do this twice, first with a large enough value of m s that the 
strange quark is heavy enough that immediately above deconfinement /x is still less 
than m s and there are still no strange quarks present. For \x = the density is 
zero; isospin and strangeness are unbroken; Lorentz symmetry is unbroken; chiral 
symmetry is broken. Above a first order transition^ at an onset chemical potential 
fi Q ~ 300 MeV, one finds nuclear matter ("nuclear" in Figure p]). Lorentz symmetry 
is broken, leaving only rotational symmetry manifest. Chiral symmetry is thought 
to be broken, although the chiral condensate (qq) is expected to be reduced from 
its vacuum value. In the nuclear matter phase, we expect an instability of the 
nucleon Fermi surfaces to lead to Cooper pairing. We assume that (as is observed 
in nuclei) the pairing is pp and nn, breaking isospin. Since there are no strange 
baryons present, U(l)s is unbroken. 

At the large value of m s we are currently describing, when /i is increased above 
jiy, we find the "2SC" phase of color-superconducting matter consisting of up and 
down quarks only, described in Refs. Jl], ^, [|, [3J]. A nucleon description is no longer 
appropriate. The light quarks pair in isosinglet channels. SU(2)y is unbroken. 
SU {2) a is unbroken. Quarks of two colors pair in a Lorentz singlet channel; quarks 
carrying the third color can form an axial vector (ferromagnetic) condensate, which 
breaks rotational invariance. The associated gap is of order a keV or much less|| 
and we neglect this condensate in Figure 1.0 The phase transition at fiy is first 



3 Discussed in Ref. 10 . 

4 In a mean- field treatment, one gluon exchange is neither attractive nor repulsive in the 
channel which leads to the ferromagnetic condensate while the instanton interaction is only 
weakly attractive. The ferromagnetic condensate is therefore fragile in the sense that this channel 
could easily be rendered repulsive by other interactions. |fj If that were the case, the third color 
quarks would presumably form Cooper pairs with nonzero orbital angular momentum, and an 
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order || f|, ||, |11|, |9], [L2| and is characterized by a competition between the chiral 



(qq) condensate and the superconducting (qq) condensate. |5|, [L2 



As the chemical potential is increased further, when [i exceeds the constituent 
strange quark mass M s (n) a strange quark Fermi surface forms, with a Fermi 
momentum far below that for the light quarks. We denote the resulting phase 
"2SC+s". Light and strange quarks do not pair with each other, because their 
respective Fermi momenta are so different (see Section |2[). The strange Fermi 
surface is presumably nevertheless unstable. The resulting ss condensate must 
be constructed from Cooper pairs which are either color 6, or have spin 1, or 
have nonzero angular momentum, or must be (sC^^s), which is symmetric in 
Dirac indices. Each of these options lead to small gaps, for different reasons: 
One-gluon exchange is repulsive in rotationally invariant color 6 channels and the 
(sC^s) channel, and although other interactions may overcome this and result 
in a net attraction, this is likely to be much weaker than the attraction in the 
dominant color 3 channels. Gaps involving Cooper pairs with nonzero J tend to be 
significantly suppressed because not all quarks at the Fermi surface participate. |3j 
Thus, although U(l)s will be broken in the "2SC+s" phase by an ss condensate in 
some channel, we expect that the resulting gaps will be very small. In our analysis 
below we therefore neglect the difference between the 2SC and 2SC+S phases. 

Finally, when the chemical potential is high enough that the Fermi momenta for 
the strange and light quarks become comparable, we cross the first order locking 
transition described in detail in the next Sections, and find the color-flavor locked 
(CFL) phase. There is an unbroken global symmetry constructed by locking the 
SU(2) V isospin rotations and an SU{2) subgroup of color. Chiral symmetry is 
once again broken. 

We now describe the sequence of phases which arise as n is increased, this 
time for a value of m s small enough that strange baryonic matter forms below 
the deconfinement density. At /i Q , one enters the nuclear matter phase, with the 
familiar nn and pp pairing at the neutron and proton Fermi surfaces breaking 
isospin. The A, £ and H densities are still zero, and strangeness is unbroken. 
At a somewhat larger chemical potential, we enter the strange baryonic matter 
phase, with Fermi surfaces for the A and S. These pair with themselves in 
spin singlets, breaking U(l)s- This phase is labelled "strange baryon" in Figure 
1. The global symmetries SU{2) L x SU{2) R and U{\)s are all broken. As /i 
rises, one finds yet another onset at which the H density becomes nonzero. This 
breaks no new symmetries, and so is not shown in the Figure. Note that kaon 
condensation fl3|1 breaks U(l)s, and SU(2)v- The only phase in the diagram in 



which this occurs is that which we have labelled the strange baryon phase. Thus, if 
kaon condensation occurs, by definition it occurs within this region of the diagram. 



even smaller gap. 
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If kaon condensation is favored, this will tend to enlarge the region of the diagram 
within which U(l)s and SU(2)v are both broken. 

We can imagine two possibilities for what happens next as \x increases further. 
(1) Deconfmement: the baryonic Fermi surface is replaced by u,d,s quark Fermi 
surfaces, which are unstable against pairing, and we enter the CFL phase, described 
above. Isospin is locked to color and SU(2) co i or+v is restored, but chiral symmetry 
remains broken. (2) No deconfmement: the Fermi momenta of all of the octet 
baryons are now similar enough that pairing between baryons with differing 
strangeness becomes possible. At this point, isospin is restored: the baryons pair 
in rotationally invariant isosinglets (pH~, nS°, AA). The interesting 

point is that scenario (1) and scenario (2) are indistinguishable. Both look like 
the "CFL" phase of the figure: U(l)s and chirality are broken, and there is an 
unbroken vector SU(2). This is the "continuity of quark and hadron matter" 
described by Schafer and Wilczek J7|. We conclude that for low enough strange 
quark mass, m s < m c s ont , there may be a region where sufficiently dense baryonic 
matter has the same symmetries as quark matter, and there need not be any 
phase transition between them. In Section 6 we use this observation to construct 
a mapping between the gaps we have calculated at the Fermi surfaces in the quark 
matter phase and gaps at the baryonic Fermi surfaces at lower densities. 

All of the qualitative features shown in Figure 1 follow from the above discussion 
of the small and large m s regimes, with one exception. In Figure 1, a transition 
between two flavor nuclear matter and three flavor quark matter (in the 2SC+S 
phase) occurs only for a range of values of m s above m c s ont . However, it may be 
that the strange baryon phase ends below m™'. One would then have a transition 
between two flavor nuclear matter and three flavor quark matter (in either the 
2SC+s phase or the CFL phase) for a range of m s values which extends both above 
and below m c s ont . Determining the extent of the strange baryon phase requires 
a detailed analysis of strange baryonic matter, including the possibility of kaon 
condensation. This is not our goal here. 

This concludes our overview of the qualitative features of the phase diagram for 
low temperature strongly interacting matter. In subsequent sections we analyze 
the 2SC and CFL phases and the unlocking transition more quantititatively, and 
use the properties of the CFL phase to make predictions for properties of baryonic 
matter in which U(l)s is broken while SU(2)y is unbroken. 

2 Model Independent Features of the Unlocking Phase 
Transition 

In this section, we give a model independent argument that the unlocking phase 
transition between the CFL and 2SC phases in Figure 1 must be first order. 
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For any m s 7^ 0, transformations in SU(3)a which involve the strange quark 
are explicitly not symmetries. The CFL and 2SC phases are distinguished by 
whether the chiral SU(2)a rotations involving only the u and d quarks are or 
are not spontaneously broken. As we will make clear in subsequent sections, the 
unlocking transition is associated with the vanishing of those diquark condensates 
which pair a strange quark with either an up or a down quark. We denote the 
resulting gaps A us for simplicity. In the absence of any A us gap, the only Cooper 
pairs are those involving pairs of light quarks, or pairs of strange quarks. The 
light quark condensate is unaffected by the strange quarks, and behaves as in a 
theory with only two flavors of quarks. Chiral symmetry is unbroken. We will 
show explicitly below that when A us 7^ 0, the interaction between the light quark 
condensates and the mixed (light and strange) condensate results in the breaking 
of two-flavor chiral symmetry SU{2)a via the locking of SU(2)l and SU{2)r flavor 
symmetries to an SU(2) subgroup of color. This color-flavor locking mechanism 
leaves a global SU (2) co i or+ y group unbroken. 

The unlocking transition is a transition between a phase with A us 7^ at 
M s < M^ nlock and a phase with A us = 0. (M s (/i) is the constituent strange quark 
mass at chemical potential /i.) The BCS mechanism guarantees superconductivity 
in the presence of an arbitrarily weak attractive interaction, and there is certainly 
an attraction between u and s quarks (with color 3) for any M s . How, then, can 
A us vanish above M" nlock ? The BCS result relies on a singularity which arises for 
pairs of fermions with zero total momentum at the Fermi surface. We see from 
Figure 2 that no pairing is possible for quarks with momenta between the u and 
s Fermi momenta, and that at most one of the quarks in a u-s Cooper pair can 
be at its respective Fermi surface. The BCS singularity therefore does not arise if 
M s 7^ 0, and a u-s condensate is not mandatory. A u-s condensate involves pairing 
of quarks with momenta within about A us of the Fermi surface, and we therefore 
expect that A us can only be nonzero if the mismatch between the up and strange 
Fermi momenta is less than or of order A us : 

vV - mm 2 - vV - m.w « MsW-mm 2 < Aug {2A) 

Here M s (fj) and M u (/i) are the constituent quark masses in the CFL phase. We 
neglect M u (fi) in the following. Equation fl2.1|) implies that arbitrarily small values 
of A us are impossible. As m s is increased from zero, A us decreases until it is 
comparable to M s (/i) 2 /2/i. At this point, smaller nonzero values of A us are not 
possible, and A us must therefore vanish discontinuously. This simple physical 
argument leads us to conclude that the unlocking phase transition at M s = j\/f" nlock 
must be first order. 

Below, we confirm by explicit calculation in a model that the unlocking phase 
transition is first order. We find that A us is of order 50 — 100 MeV on the CFL 
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Figure 2: How the strange quark mass disrupts a u-s condensate. The strange 
quark (upper curve) and light quark (straight line) dispersion relations are shown, 
with their Fermi seas filled up to the Fermi energy Ep. The horizontal axis is the 
magnitude of the spatial momentum; pairing occurs between particles (or holes) 
with the same p and opposite p. For p < p s F , hole-hole pairing (s-u) is possible 
(two examples are shown). For p > p F , particle-particle pairing (s-u) is possible 
(one example is shown). Between the Fermi momenta for the s and u quarks, no 
such pairing is possible. 

side of the transition if the coupling is calibrated to give a reasonable magnitude 
for the chiral condensate in vacuum. 

3 Superconducting Condensates in a Model 

We study the physics of the quark matter phases in Figure 1 in a toy model in 
which we replace the full interactions between quarks by a four-fermion interaction 
with the quantum numbers of single-gluon exchange, 

interaction = G j d\ (q\ A -fq) (qX^q) , (3.1) 

and work in a mean-field approximation. All four-fermion interactions involving 
fermions at the Fermi surface are equally relevant, in the renormalization group 
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sense, so why use only the one with the quantum numbers of one-gluon exchange? 
Renormalization group analyses fll4], |15|| show that other interactions are important, 



but confirm that in QCD with two and three massless quarks the most attractive 
channels for condensation are those generated by ( |3.1| ). There is one important 
caveat. The single-gluon-exchange interaction is symmetric under U(1)a, and so it 
sees no distinction between condensates of the form (qCq) and (qC^q). However, 
once instantons are included the Lorentz scalar (qCj^q) is favored, 0, |] so we 
neglect the other form. 

One of the central qualitative lessons of past work is that superconducting 
gaps of order 100 MeV are obtained in the quark matter phase independent of the 
details of the interaction which is used, as long as the strength of the interaction is 
chosen in a way which is roughly consistent with what we know about the vacuum 
chiral condensate. Superconducting gaps of this magnitude can be obtained 
using a number of different treatments based upon single-gluon exchange. 0, [16], || 



Refs. ||, H [12] find gaps of this magnitude by making approximations of varying 
sophistication in which the interaction between quarks is modeled by that induced 
by instantons. Further confirmation of this lesson comes from recent work using 
the NJL model. flTTI 



As in any model which uses a four-fermion interaction, in addition to choosing 
the quantum numbers of the interaction we must introduce an ad hoc form factor. 
The explicit lesson of Refs. [0, ||, |17|, implicit in other work, is that results are 
rather insensitive to the choice of form factor, again as long as G is suitably 
normalized. We therefore use a simple sharp cutoff in momentum integrals at 
a spatial momentum A = 800 MeV. Our results should therefore not be extended 
above /i ~ A. The form factor (which we describe by the single parameter A) is a 
crude representation of physics not described by the model, and A should not be 
taken to infinity. For A = 0.8 GeV, setting G = 7.5 GeV~ 2 in (|3.1| ) results in a 



vacuum constituent mass for the light quarks of 400 MeV, and we use these values 
of A and G throughout. We have however checked that if we instead take A = 1. 
GeV, and change G accordingly, the superconducting gaps and the location of the 
unlocking phase transition are not significantly affected. 
An analysis of the phase diagram requires condensates 

(tfCTM?), (^7574^), (fig?), (3.2) 
leading to gap parameters 

Ag, "g, <• (3-3) 

with corresponding quantum numbers. Each gap matrix is a symmetric 9x9 
matrix describing the color (Greek indices) and flavor (Roman indices) structure. 
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In |Appendix A] we derive the gap equation for these condensates in full 



generality, assuming nothing about their color-flavor structure. However, it is 
very difficult to solve the general gap equation, so at this point we make several 
simplifying assumptions to obtain easily soluble gap equations. These are: 

1. We fix = M s 5 l 3 5]5P. For more details see below. 

2. We use the simplest form for A"- which allows an interpolation between 
the color-flavor locking favored by single-gluon exchange at m s = and the 



"2SC" phase favored at m s — > oo (see ( |3.4j) ). This ansatz, which requires 
five independent superconducting gap parameters, leads to consistent gap 
equations in the presence of the one-gluon exchange interaction. We describe 
this ansatz in detail in Section 4. 

We neglect . This condensate pairs left-handed and right-handed quarks, 
and so breaks chiral symmetry. It can be shown to vanish in the absence 
of quark masses. [T^] In [Appendix A| we show that it must be nonzero in the 



presence of a nonzero (qCj^q) condensate and nonzero quark masses. In 
the 2SC phase, the (qC^q) condensate involves only the massless quarks 
and (qC'js'y^q) therefore vanishes and chiral symmetry remains unbroken. 
In the CFL phase, however, the (qC'y^q) condensate involves the massive 
strange quarks and itself breaks chiral symmetry by color-flavor locking. No 
symmetry argument precludes k^O, and in [Appendix A| we have confirmed 



by direct calculation at one value of /i that the gaps are nonzero in the 
CFL phase. However, we find that these gaps are much smaller than the 
corresponding gaps generated by the (qC^q) condensate. Including k in the 
gap equations for A makes finding solutions prohibitively slow. Therefore, 
once we have convinced ourselves that these condensates are present but 
small, we neglect them. 

We now give a more detailed discussion of our assumptions about the quark 
mass matrix 0. Since all our calculations are in the quark matter regime (yU > fiy 
in Figure |l|), we will assume that /i is high enough that we can neglect chiral 
condensates for the u and d quarks^ but because the current quark mass m s 
is nonzero there may be a nonzero (ss) chiral condensate. In this case the quark 
mass matrix is = diag(0, 0, M a ), where M s is the /z-dependent constituent strange 
quark mass, satisfying a gap equation of its own. At sufficiently high densities, it 

5 Because chiral symmetry is broken in the CFL phase, there is no reason for the ordinary 
(uu) and (dd) chiral condensate to vanish. Indeed, Ref. || demonstrates that the (uu) and (dd) 
condensates must be nonvanishing. However, an explicit calculation |18|] shows them to be small, 



and we neglect them. Eq. (2.1) shows that nonzero (uu) and (dd) condensates allow color-flavor 



locking to persist to higher strange quark masses. 
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is given by the current mass m s , of order 100 MeV. At lower densities, it receives 
an additional contribution from (ss). In the interests of simplicity, however, we 
will not solve the M s gap equation simultaneously with the superconductivity gap 
equations. Instead, we treat M s as a parameter in the superconductivity gap 
equations, and do not determine what value of M s corresponds to given values 
of m s and \i. For simplicity, we choose the same value of the parameter M s for 
strange quarks of all three colors. This would not be the case in a treatment in 
which the M s gap equations are solved simultaneously with the superconductivity 
gap equations. [|l| 



With our simplifying assumptions, we are left needing only a gap equation for 

iterp< 
0) is 



the quark-quark condensate A"^. The simplest ansatz that interpolates between 



the two flavor case (m s 



/ b + e 
b 

c 



A af3 



b c 
b + e c 
c d 



oo) and the three flavor case (m s 



f 



f 



\ 



f 

f J 



(3.4) 



basis vectors: 

(ar,i) = (1,1), (2,2), (3,3), (1,2), (2,1), (1,3), (3,1), (2,3), (3,2) 
= (r,u), (g,d), (b,s), (r,d), (g,u), (r,s), (b,u), (g,s), (b,d) 

where the color indices are a,/3 and the flavor indices are The strange quark 
is i = 3. The rows are labelled by (a,i) and the columns by (/3,j). 



description 



condensate 



symmetry 



2SC: 


2-flavor 

sup er conduct i vity 


c = 
b = 


d = f = 0, 
—e 


SU(2) L x 


SU(2) R 


CFL 


color-flavor locking 






SU (2) color 


i-L+R 


CFL 


color-flavor locking 
with m s = 


c = 
d = 


b,f = e, 
b + e 


SU (3) co i or 


fL+-R 



Table 1: Symmetries of the condensate ansatz (|3.4j) in various regimes 



The properties of the ansatz are summarized in Table [|. In its general form, 
this condensate locks color and flavor. This is because of the condensates c and 
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/, referred to collectively as A us above, that combine a strange quark with a light 
one. It is straightforward to confirm by direct calculation that if either c or / or 
b + e is nonzero, then the matrix A?- of (|3.4|) is not invariant under separate flavor 
or color rotations but is left invariant by simultaneous rotations of SU (2)y and the 
SU(2) subgroup of color corresponding to the colors 1 and 2. Thus, color-flavor 
locking occurs whenever one or more of c, /, or b + e is nonzero. We discuss b + e 
below. 

Although the standard electromagnetic symmetry is broken in the CFL phase, 
as are all the color gauge symmetries, there is a combination of electromagnetic and 
color symmetry that is preserved. Consider the gauged U(l) under which the 
charge Q' of each quark is the sum of its electromagnetic charge (2/3,-1/3,-1/3) 
(depending on the flavor of the quark) and its color hypercharge (—2/3, 1/3, 1/3) 
(depending on the color of the quark). It is easy to confirm that the sum of the 
Q' charges of each pair of quarks corresponding to a nonzero entry in ( |3.4j ) is zero. 
This modified electromagnetism is therefore not broken by the condensate. 

At M s = 0, one has color-flavor locking: c, / and b + e are all nonzero. In 
addition, c = b, f = e, and d = b+e, and the matrix is invariant under simultaneous 
rotations of SU(3)y and S77(3) co i r- For M s nonzero but sufficiently small, c, / 
and b + e remain nonzero but are no longer equal to b, e and d. Color and flavor 
are locked, and the matrix is invariant under simultaneous rotations of SU(2)v 
and S77(2) co i or . As described in Section 2, once M s becomes large enough, c and / 
both vanish (and so does b + e) and the symmetry group enlarges, unlocking color 
and flavor, and restoring chiral symmetry (see Table |l|). 

If A were nonzero only in color 3 channels, we would have b = — e, d = and 
c = — /. Single-gluon-exchange is repulsive in the color 6 channels. However, in 
the CFL phase where c and / are nonzero, the only consistent solutions to the gap 
equations have b ^ — e, c ^ —f, and d ^ 0. This means that the single-gluon- 
exchange interaction requires a small color 6 admixture along with the favored color 
3 condensate, even though it prohibits color 6 condensates alone. This feature 
occurs both for M s = 0[§] and for < M s < M s unlock . Note that not all color 6 
terms arise in ( |3.4j) . For example, only three of the nine diagonal elements of Q3.4Q 
are nonzero. These are the only three diagonal elements which can be nonzero 
without breaking the SU(2) co i OT+v symmetry, and upsetting color-flavor locking. 
Furthermore, if the two upper-most diagonal elements were not equal to b + e, the 
SU(2) color+v symmetry would also be broken. We have therefore learned that the 
color 6 terms which are induced in the presence of a color 3 condensate are those 
which do not change the global symmetry of the color 3 condensate. Those color 
6 terms which are "allowed" in this sense are generated by the interaction. 

The corresponding analysis of the 2SC phase, where only two flavors participate 
in the condensate is similar in logic but leads to different conclusions. In the 2SC 
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phase, where c and / are zero, we do in fact find b = — e and d = 0, namely a pure 
color 3 condensate. In this case, if b + e were nonzero, one would have color-flavor 
locking and broken chiral symmetry. In this case, then, the pure color 3 condensate 
with b + e = is "protected" by the fact that any admixture of color 6 condensate 
like b + e would break a global symmetry which the color 3 condensate does not 
break. In the absence of c and /, a nonzero d condensate ((ss)) is still allowed, 
but there is no way for the color 3 condensate to induce it. Therefore the simple 
argument that one-gluon exchange is repulsive in this channel holds, and in our 
model we find d = in the 2SC phase even where the strange quark density is 
nonzero. As discussed above, a more complete analysis would include a small (ss) 
condensate in some channel, resulting in a distinction between the 2SC and 2SC+S 
phases which our analysis does not see. 



4 The Gap Equation . . . 

Our task now is to obtain the gap equation for A, and verify the behavior described 
in Sections 1 and 2. In particular, we want to show that the unlocking phase 
transition is first-order, and calculate M" nlock (/i) and compare to the criterion 
( |2.1| ) which we derived on model-independent grounds. 

In |Appcndix A we have derived the mean field gap equation and given its 



general solution. It takes the form (see Eq. ( |A.5|) ) 

A = J^i Jd 4 q\ T Al,Pi(^M s ,A,q)\ A ^, (4.1) 



where the function Pi is given by (TADD and flAJD , and we are assuming, as 
described in Section [3] that 0^- = M s S l 3 SjS^ and that k = 0. Note that strictly 
speaking the gap equations do not close for k — 0. The k gap equation is just 
like ( |4.1| ) but with — P 4 instead of P\ on the right hand side, and for M s > 0, 
P4 is nonzero even when k = 0. This means that it is inconsistent to set k = 0. 
However, as discussed in the Appendix, k turns out to be small and we use ( [4.1|) 
with k = in Pi. 

In principle we could just evaluate Pi(/i, M s , A, q) by substituting (|3.4| ) into 
(|A.7|) , but that would lead to a very complicated expression. To turn (|4.1| ) into a 
tractable set of gap equations, we transform to a slightly different basis from ( p.4| ). 

In the new basis, the first 2 vectors are (a, i) = l/y/2,((l, 1) ± (2, 2)V Now the 

top left 3x3 block of (|J) looks like 



(4.2) 









r 

an 


an 


H 


\ a 12 


a-22 J 
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A is now block-diagonal, consisting of 1 x 1 and 2x2 blocks. We find that 
P\(fi, M s , A, q) takes a corresponding block-diagonal form, 



Pi(/i,M„A,g) 



/ E(e,q) 

A n (a,q) A 12 (a,q) 
A 12 (a,q) A 22 (a,q) 



E(e,q) 



E(e,q) 



F(f,q) 



F(f,q) 



F(f,q) 



F(f,q) 



(4.3) 



where 



F(f,q) = P 1 ( t i,M s ,f,q) 



fw 



E(e,q) = Pi(jU, 0,e,q) 
A{a,q) = PtifJ,, M s , a, q) 



w 2 - (4/i 2 - M 2 ) q 2 - (/i - iq ) 2 M 2 
w = f 2 + fi 2 + q 2 + q 2 ) 

as above, with / — > e and M s = 



(4.4) 



Note that the function Pi simplifies considerably when its argument is a symmetric 
off-diagonal 2x2 matrix, as in the e and / blocks. The resulting expressions for 
E and F are given on the right-hand side of ( (4.4| ) . The only difficult thing in (O) 
is evaluating Pi for the 2x2 matrix a, in which case Eqns. ( |A.6| ) and (|A.7| ) must 
be used. 

The only nonzero matrix elements in Pi, on the right hand side of the gap 
equation, are those which are nonzero in our ansatz for A, on the left hand side of 
the gap equation. This confirms that the only color 6 terms which are induced by 
the color 3 condensate are those which we have included in our ansatz, in agreement 
with the symmetry arguments given above. Although more complicated ansatze 
may be worth considering in future work, this ansatz is the minimal one which 
suffices. 

Taking into account the A and 7 matrices from the gluon vertex, we end up 
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with the 5 gap equations 



e 


= 4G/(2vr 


) 4 /A 


(-A u (a,g) - §£(e,g)) 


/ 


= 4G/(2vr 


) 4 /A 


(V2A 12 (a, q) - §F(/, g)) 


a n 


= 4G/(2tt 




(iA n (a,g) + 3E(e,g)) 


an 


= 4G/(2tt 


) 4 / A 


(-|A 12 (a,g)+2 v / 2F(/,g)) 


a-22 


= 4G/(2tt 


) 4 /A 


§A 22 (a,g) 



(4.5) 



To obtain the gap parameters as a function of fi and M s , we solve the 5 
simultaneous equations ( |4.5| ) numerically, using Mathematica. The right-hand sides 
are evaluated by numerical integration of the integrands given in (|4.4| ). For A(a, q), 
we must evaluate M s , a, q), which corresponds to (|A.7|) with A being the 2x2 
matrix a, k = 0, and = diag(0,M s ). The gap parameters b,c,d,e,f determine 
the poles of the propagator Pi and hence determine the quasiparticle dispersion 
relations and the gaps at the Fermi surface. 

5 ... and Its Solutions 



In this section, we present solutions to the gap equations ( |4.5| ). We demonstrate 
explicitly that the unlocking phase transition is first order in our model, as we have 
argued on general grounds in Section 2. Furthermore, we confirm that the simple 
criterion (|2.1|) is a good guide. We determine the value of M™ nt , the highest strange 
quark constituent mass at which, in the model, no 2SC phase intrudes. We use 
our results to draw tentative conclusions for the behavior of strongly interacting 
matter function of density in nature. 

We plot the solutions b, c, d, e, / to the gap equations along two different lines of 
constant M s and one line of constant fi, all shown in Figure |^. Recall that b = —e 
and c = d = / = 0in the 2SC phase whereas all the gaps are nonzero in the CFL 
phase. For M s = 0, b = c and e = /. If the condensate were purely color 3 in the 
CFL phase, one would have b = — e, c = — / and d = 0. 

Figures 4 and 5 show the solutions to the gap equations as a function of fi, 
for M s = 150 MeV and M s = 350 MeV. We see that for M s = 150 MeV, the 
color-flavor-locked phase continues down to fi = 400 MeV, below which we expect 
a baryonic description to become appropriate. In contrast, in Figure 5 we see 
that for M s = 350 MeV, the CFL phase only exists for ^ > 530 MeV. At lower 
densities, below a first order unlocking phase transition, the 2SC phase is favored. f] 
This transition occurs at high enough densities that a quark matter description is 



3 In the calculations for our plots we assume M s is the same on either side of the phase 
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400 MeV 



Figure 3: The part of the phase diagram in which we have done calculations using 
our model. The thin solid lines show the slices along which we have solved the gap 
equations. Note that a horizontal line in Figure 1 at some m s corresponds to a 
curve in this figure which approaches the horizontal line M s = m s from above as 
li — > oo. 



justified, and our model can therefore be used to describe it. At lower densities still, 
there is another first order phase transition to a baryonic phase. 0, f|, [j[ [□], |9|, [12 



In order to estimate M^ ont , in Figure 6 we plot the solutions to the gap equations 
as a function of M s with \i fixed at 400 MeV. This is the lowest fi at which we 
expect a quark matter description to be valid. We find that the CFL phase exists 
for M s below a first order transition at M^ nlock = 235 MeV. Our model therefore 
suggests that M s cont is of order 250 MeV. 

The question, then, is whether Figure 4 or Figure 5 is closer to the behavior of 
strongly interacting matter as a function of density in nature. If Figure 5 shows the 
correct qualitative behavior, we would expect two first order phase transitions as 
the density is increased above nuclear density. After the first transition, one would 
have quark matter in the 2SC phase. In this quark matter phase, the constituent 



transition. In reality M s will be somewhat smaller in the 2SC phase because there is no chiral 
symmetry breaking there, but since the 2SC condensates do not involve the s quark this is of no 
consequence. 
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Figure 4: Solutions to the gap equations as a function of /i, at M s = 150 MeV. 
At this value of M s , the CFL phase persists down to /x = 400 MeV, and below. 
The two solid lines are b and — e, the dashed lines are c and — /. At large /i, the 
gaps approach the M s = pattern. 

strange quark mass must be greater than M^ ont ~ 250 MeV. Then, above the 
next transition, a color-flavor-locked quark matter phase is obtained. If Figure 4 is 
the better gulide to nature, as we surmise, then there is no window of /i in which 
the 2SC phase intrudes, and nature may choose a continuous transition between 
baryonic and quark matter with the symmetries of the CFL phase. 

We have analyzed the unlocking phase transition quantitatively in our model, 
and can now confirm the predictions of the model-independent but qualitative 
arguments of Section 2. We first demonstrate that the transition is first order, and 
justify the location of the discontinuities shown in Figures 5 and 6. To find the 
first order transition in our model, we need the effective potential Q. We use the 
fact that the gap equations are equivalent to the requirement that the derivative 
of fl with respect to the gap parameters vanishes. The set of gap equations ( fL5|) 
corresponds to five differential equations dfl/dl | m i n = where / = (b, c, d, e, /) at 
the potential minimum, so we have analytic expressions for dQ/dl for any value 
of the gap parameters I. At a first order phase transition, Q should have two 
degenerate minima, one corresponding to the CFL phase and the other to the 
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Figure 5: Solutions to the gap equations as a function of //, at M s = 350 MeV. 

At this large value of M s , the CFL phase exists only for \x > 530 MeV. That is, 
M uniock = 350 MeV for ^ = 530 MeV _ 

2SC phase. We have found two minima at M s = 235 MeV, /x = 400 MeV. (See 
Figure 6.) We define a straight line l(x) in the five- dimensional space of gap 
parameters which goes from the 2SC minimum Z(0) = (b = 100, c = 0, d — 0, e = 
-100,/ = 0) MeV to the CFL minimum 1(1) = (91,56,4,-88,-52) MeV. We 
evaluate (dQ/dl ) along this line, and then obtain the potential itself by integrating 
(dQ/dl)(dl/dx) with respect to x. The result is shown in Figure [7|. The phase 
transition is first order. 

The model-independent arguments of Section 2 lead to the criterion (|2.1j) for 
the location of the unlocking phase transition. We now compare this prediction to 
the quantitative results we have obtained in our model. For \i = 400 MeV, at the 
critical strange quark mass M^ nlock = 235 MeV (see Figure 6) the gaps c and — /, 
which vanish in the 2SC phase, are within a few MeV of each other in the CFL 
phase, with c ~ — / ~ 55MeV. Taking A us m 55 MeV, we find 

( M unlock)2 w 2.5/iA us . (5.1) 
The unlocking phase transition in Figure 5 occurs at M s = 350 MeV, \x = 530 MeV 
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Figure 6: Solutions to the gap equations as a function of M s , at fi = 400 MeV, 
along the vertical line in Figure 3. When the strange quark becomes massive 
enough, s-u condensates become too expensive, and there is a first-order phase 
transition from CFL to 2SC. 

and has c ~ — / ~ 85 MeV; this yields a relation as in (|5.1| ), with 2.5 replaced 
by 2.7. We conclude that fl2.1|), derived by simple physical arguments based on 
comparing the mismatch between the up and strange Fermi momenta with the A us 
gap, is a very good guide to the location of the unlocking phase transition. 

6 Quark-hadron continuity 

As has been emphasized in Section |], for low enough m s the CFL phase may 
consist of hadronic matter at low fi, and quark matter at high fi. This raises the 
exciting possibility |7j] that properties of sufficiently dense hadronic matter could 
be found by extrapolation from the quark matter regime where models like the one 
considered in this paper can be used as a guide at moderate densities, and where 
the QCD gauge coupling becomes small at very high densities. 

The most straightforward application of this idea is to relate the quark/gluon 
description of the spectrum to the hadron description of the spectrum in the CFL 
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Figure 7: The effective potential Q in units of 10~ 5 GeV 4 in the vicinity of 
the unlocking phase transition, for fi = 400 MeV, M s = 235 MeV. The plot 
corresponds to a slice of the potential in the 5 dimensional parameter space along 
a line which linearly interpolates between the minima in the 2SC phase (x — 0) 
and the CFL phase (x = 1). The phase transition is clearly first order. 

phase. M As fi is decreased from the regime in which a quark/gluon description 
is convenient to one in which a baryonic description is convenient, there is no 
change in symmetry so there need be no transition: the spectrum of the theory 
may change continuously. Under this mapping, the massive gluons in the CFL 
phase map to the octet of vector bosons;Q the Goldstone bosons associated with 
chiral symmetry breaking in the CFL phase map to the pions; and the quarks map 
onto baryons. Pairing occurs at the Fermi surfaces, and we therefore expect the 
gap parameters in the various quark channels, calculated in Section ^j, to map to 
the gap parameters due to baryon pairing. 

In Table |2| we show how this works for the fermionic states in 2+1 flavor QCD. 
There are nine states in the quark matter phase. We show how they transform 
under the unbroken "isospin" of S , f/(2) color+ y and their charges under the unbroken 
"rotated electromagnetism" generated by Q', as described in Section 4. Table § also 
shows the baryon octet, and their transformation properties under the symmetries 
of isospin and electromagnetism that are unbroken in sufficiently dense hadronic 
matter. Clearly there is a correspondence between the two sets of particles.^ As 

7 The singlet vector boson in the hadronic phase does not correspond to a massive gluon in 
the CFL phase. This has been discussed in Ref. (7). 

8 The one exception is the final isosinglet. In the jj, — > 00 limit, where the full 3-flavor symmetry 
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Quark 



SU (2) co i or+v Q' gap 



Hadron SU(2)y Q gap 



bu 
bd 

gs 
rs 



+ 1 



-1 



/ 




+ 1 





-1 



A B 



/ ru — gd \ 

gu 

\ rd j 





+1 
-1 





+1 
-1 



A B 

^3 



ru + gd + £_&s 







A 







ru + gd — £ + bs 



A_ 



I ( 2b + e + d ± y/(2b + e - d) 2 + 8c 2 
+ v 7 ^ + e - d) 2 + 8c 2 



Table 2: Comparison of states and gap parameters in high density quark and 
hadronic matter. 



fi increases, the spectrum described in Table 2 may evolve continuously even as 
the language used to describe it changes from baryons, SU{2)y and Q to quarks, 
SU(2) co i or+v and Q'. 

If the spectrum changes continuously, then in particular so must the gaps. As 
discussed above, and displayed explicitly in (|3.4|), the quarks pair into rotationally 
invariant, Q'-neutral, SU(2) co \ or+ v singlets. The two doublets of Table || pair with 
each other, with gap parameter /. The triplet pairs with itself, with gap parameter 
e. Finally, the two singlets pair with themselves. 

When we map the quark states onto baryonic states, we can predict that 
the baryonic pairing scheme that will occur is the one conjectured in Sect. [I] for 
sufficiently dense baryonic matter: 



( P E-),(E-p),(nE°),(E n) 
(£+£-), (£-£+), 
(AA) 



B 



4 quasiparticles, with gap parameter A 
3 quasiparticles, with gap parameter A 3 
1 quasiparticle, with gap parameter Af 



B 



(6.1) 



is restored, it becomes an SU (3) singlet, so it is not expected to map to any member of the baryon 
octet. We discuss this further below. The gap A + in this channel is twice as large as the others 
(it corresponds to Ai in Ref. B). 
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The baryon pairs are rotationally-invariant, Q-neutral, SU(2)v singlets. It 
seems reasonable to conclude that as /i is increased the baryonic gap parameters 
(Af , Af , Af ) may evolve continuously to become the quark matter gap parameters 
(/, e,A_), which we have calculated in this paper. Assuming continuity, the 
magnitude of the gaps will change as the density is increased but if their ratios 
change less we can use our results for the gap parameters in the CFL phase at 
H = 400 MeV and M s = 150 MeV to suggest baryonic gap parameters with ratios 

Af : Af : Af ~ 1.06 : 1.26 : 1 (6.2) 

in matter which is sufficiently dense but still conveniently described as baryonic. 

As mentioned above, the ninth quark corresponds to a singlet baryon which 
is very heavy, for reasons which have nothing to do with our considerations. We 
therefore expect that as fi is increased and a quark matter description takes over 
from a baryonic description, the density of the singlet quark/baryon becomes 
nonzero at some /i, and begins to increase. Deep into the quark matter phase, 
there is a gap A + for this ninth quark, but this does not correspond to any gap 
deep in the baryonic phase, because the density of the corresponding baryon is 
zero. There is no change in symmetry at the \i at which the density of the ninth 
quark/baryon becomes nonzero, just as there was no change in symmetry at the 
lower fi (within the "strange baryon" phase of Figure 1) at which the S density 
became nonzero. Both these onsets are smooth transitions at arbitrarily small 
but nonzero temperature. We do not expect the onset of a nonzero density for 
the ninth quark/baryon to upset the continuity between the baryonic and quark 
matter phases, but symmetry arguments can only demonstrate the possibility of 
continuity; they cannot prove that there is no transition. 

The analysis leading to ( |6.2|) provides the first example of the use of the 
hypothesis of quark-hadron continuity in the color-flavor locked phase to map a 
quark matter calculation onto quantitative properties of baryonic matter. 



7 Color-Flavor Locking at Asymptotic Densities 

Son P| has studied color superconductivity at much higher densities than those 
we can treat with our model, using a resummed gluon propagator rather than a 
point-like four-quark interaction, and has determined the leading behavior of the 
gap in the large /z, small g(/i) limit. Here, g(q) is the QCD gauge coupling at 
momentum transfer q. The resulting expression, 

1 ( 3tt 2 1 \ 

describes the /i-dependence of the gap(s), but not their absolute normalization. A 
could be any of our gaps. The distinction between b, c, d, e and / is in the prefactor 
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Figure 8: The upper curve shows Son's result for the superconducting gap A as 
a function of log 10 /i, for \i from 0.8 GeV to 10 6 GeV. The vertical scale has been 
normalized so that A = 0.1 GeV at fi — 0.8 GeV. We have taken g(fi) from 
the two-loop beta function for three flavor QCD with Aqcd = 200 MeV. Color- 
flavor locking occurs whenever A > M^/2/i. The lower curve is M^/2/i, taking 
M s = 150 MeV. We conclude that QCD at very high densities is in the CFL phase. 

which should appear in ( |7.1| ), but which has so far not been calculated. We take 
our model estimates as a guide at moderate densities, around \i ~ 400 — 800 MeV, 
and use them to normalize the calculation of Ref. M by fixing the prefactor in ( |7.1| ) 
so that A = 100 MeV at /i = 800 MeV. We show the result in Figure §. Note that 
A is plotted versus log/i; it changes very slowly. It decreases by about a factor 
of three as /x is increased to around 100 GeV (!) and then begins to rise without 
bound at even higher densities. We conclude that independent of any details (like 
the precise value of M s , for example) at asymptotically high densities A is far 
above M s 2 /2/z. For any finite value of the strange quark mass m s , quark matter is 
in the color-flavor locked phase, with broken chiral symmetry, at arbitrarily high 
densities where the gauge coupling becomes small. 

8 Conclusions 

We have discussed a conjectured phase diagram for 2 + 1 flavor QCD as a function 
of /i, the chemical potential for quark number, and m s , the strange quark mass. 
We have worked at zero temperature and ignored electromagnetism and the u and 
d quark masses throughout. Our quantitative analysis is restricted to the high 
density quark matter phases, where we use a model with an effective four-fermion 
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interaction which has the index structure of single gluon exchange. 

Our analysis shows that for any finite m s the quarks will always pair in a color- 
flavor-locked (CFL) state at arbitrarily large chemical potential, breaking chiral 
symmetry. What happens at lower \x depends on the strange quark mass. If m s is 
greater than a critical value m c s ont , then as \i is reduced, there will be an unlocking 
phase transition to a color superconductor phase analogous to that in two flavor 
QCD , with restoration of chiral symmetry. We give model-independent arguments 
that the unlocking transition between the CFL and 2SC phases is controlled by 
the mismatch between the light and strange quark Fermi momenta and is therfore 
first order. We confirm this quantitatively in our model. If \x is reduced further, 
there will presumably be a phase transition to nuclear matter, and chiral symmetry 
will be broken again. However, if m s < m c ° nt then the 2SC state never occurs. 
The quark matter stays in the color-flavor locked phase all the way down until the 
transition to hadronic matter. Chiral symmetry is never restored at any [i. In this 
case there may be continuity between the CFL quark matter and sufficiently dense 
hadronic matter. Assuming such continuity, we use our calculation of quark gaps 
to make predictions about the gaps in hadronic matter ( |6.2|) . 

At arbitrarily high densities, where the QCD gauge coupling is small, quark 
matter is always in the CFL phase with broken chiral symmetry. This is true 
independent of whether the "transition" to quark matter is continuous (as may 
occur for small m s , including, we surmise, realistic m s ) or whether, as for larger 
m s , there are two first order transitions, from nuclear matter to the 2SC phase, 
and then to the CFL phase. 

There are many directions in which this work can be developed. We have 
worked at zero temperature, so a natural extension would be to study the 
effects of finite temperature. The phase diagram of two-flavor QCD as a 
function of baryon density, temperature and quark mass has been explored in 
Ref. ||. A nonperturbative approach beyond the mean field approximation that 
we have employed can be performed along the lines of Ref. |2D| using the exact 
renormalization group, or by doing a lattice calculation ||21|| . 

Within our model, one could study more exotic channels that we have ignored, 
such as those with S = 1 and/or L — 1, or channels that would lead to pairing of 
the strange quarks in the 2SC+S phase. There are also improvements that could 
be made to the model, such as including four- and six-fermion interactions induced 
by the three-flavor instanton vertex, or using sum rules to include nonperturbative 



gluons p2[ , or including perturbative gluons ||. It would also be desirable to 
include the effects of electromagnetism, and of different chemical potentials for the 
different flavors. 

Finally, it is of great importance to investigate the consequences of our findings 
for the phenomenology of neutron/quark stars, which are the only naturally 
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occuring example of cold matter at the densities we have studied. We are confident 
that, with the basic symmetry properties of the phase diagram now at hand, a 
whole new phenomenology waits to be uncovered as the role of the strange quark 
in dense matter is fully elucidated. 
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Appendix A The general form of the gap equation 

We calculate the gap equations for a quark-quark condensate with completely 
general color-flavor structure, using a four-fermion interaction vertex with the 
index structure of single-gluon exchange. We will allow condensation in the 
channels A^ 8 = (qfC^qj) and nff = (q^C^j^qj), with arbitrary color-flavor 
structure. We will also allow the quark mass matrix <fi?jr to have arbitrary color- 
flavor structure. We use the Euclidean conventions of Ref. || 



The mean field gap equation is 

G 



(27T) 



Jd'qV^M- 1 ^, (A.i; 



Since we want to study quark-quark condensation, we use the Nambu-Gorkov basis, 
in which the quark-gluon vertex has the standard index structure 

V A , = ( , „ ) , (A.2) 



and the fermion matrix M is 

M(q) 



(A.3) 



(g ± if^h + Qil 1 



C is the Dirac charge conjugation matrix. In this appendix we allow A, k, to 
have arbitrary color-flavor structure, so we treat them as general non-commuting 
matrices in 

C>y 5 (A + K<y 4 ) -i(f) 

up (A + kj^C 



(A.4) 



25 



We will fix the quark masses 0, and solve for the superconducting gaps A and k. 
We can just take the top left component of (|A.1|) , which gives us 

C 75 (A + K l4 ) = JLJ d*q\ T Al l X- 1 \ Al », (A.5) 
where X' 1 = (M~ 1 ) 22 , and X can be written 
X = 75 (AV + X 47 4 + Xvqrf + Xsqia^C 

X 1 = (q 2 + fi 2 )R + (J)R(f) + i(qo-ilj)A- 1 KR(j)-i(qo + ifi)(j)A- 1 KR + A 
X 4 = i(q -in)R(j) - i(q + ifj)<f>R + (q 2 - q 2 + (i 2 )A~ 1 kR + (pA^nR^ - « 
ATy = 2q A- 1 kR + i{R(j> - <f>R) 
X s = 2fiR + A- 1 nR(p + (pA^nR 
R = (A-kA^k)- 1 . 



(A.6) 



X can now be inverted, 



- 1 = -C^PJ + P^ + Pvqa' + Ps^) 

D i = -(U7 1 u 1 -Tr 1 Ti)- 1 ur 1 sr 1 



X- 1 = -C(p 1 I + P 4l i + P v q il i + P s q i a^) l5 

Pi 

Pi 
Lh 

U A = -S 2 -\-X 1 + q<X s X^X v ) + S7\iq £ X v X^X s -X A ) (A 7) 
„i, ■ -« - - - - 1 - - - - 1 - .----1--, 

J 1 
T 4 

^ = <f 2 X 5 + iq 2 X v Xr 1 X 4 

5 2 = q 2 (iX v + X s X^ l X 4 ) 

53 = X\ — x 4 x 1 x 4 

We have ignored Py and P$ because these multiply odd powers of momentum, and 
so give no contribution to the integral on the RHS of the gap equation ( |A.5| ). 




Note that in general P 4 7^ 0, so it is necessary to include k in the ansatz. 
Without it the gap equations do not close. That is, if one sets k to zero, one finds 
that if M s and A are both nonzero, one obtains P 4 7^ and concludes that setting 
k = is inconsistent. However, k turns out to be small, and we neglect it. To 
check that k is small, we solved the coupled gap equations for the k and / gaps, in 
one of the 2x2 blocks of the ansatz (|3.4| ), using a modified gluon vertex so that 
this block did not mix with the others. We calculated k and / as a function of 
M s , at n — 400 MeV. The values of / were very close to those obtained for the 
full gluon vertex (see Figure H), and k was zero at M s = 0, and always less than 
6 MeV. 
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